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Математичне моделювання в систем! автоматичного керування 


процесом пневмогдродинамично! ди на шдземн! масиви 


Рассматривается задача построения системы автоматического управления технологическим процессом 
пневмо- и гидродинамического воздействия на газонасыщенный подземный угольно-породный массив 
как способа направленного изменения его состояния. Для описания процесса предлагается матема- 
тическая модель, в основу которой положена система нелинейных уравнений с краевыми условиями, 
определяемыми технологической схемой. 
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Розглядаеться завдання побудови системи автоматичного керування технолог!чним процесом пневмо- та 
пдродинам!чно! дй на газонасичений шдземний вупльно-пордний масив як засобу спрямовано! змни 
його стану. Для опису процесу пропонуеться математична модель, в основу яко! покладено систему 
нелинйних р1внянь 3 крайовими умовами, що визначаються технолог!чною схемою. 
Ключов! слова: процес, математична модель, керування, масив, система. 


Введение 


Актуальность задачи. При разработке месторождений полезных ископаемых, 
в частности, газоносных пластов, в процессе эксплуатации месторождения возникает 
проблема повышения интенсивности извлечения полезного продукта. Одним ИЗ ОСНОВ- 
ных способов решения проблемы является целенаправленное изменение свойств мас- 
сива путем внешнего воздействия, в частности, пневмогидродинамической обработки. 
Качественная реализация воздействия обеспечивается с помощью автоматического уп- 
равления его параметрами, при этом состояние процесса представляется с помощью 
математической модели. При проектировании системы управления возникает необхо- 
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димость использования элементов искусственного интеллекта, в частности, элементов 
нечеткой логики. В этой связи тема работы является актуальной. 

Цель работы — обоснование структуры математической модели и алгоритма 
управления процессом пневмогидродинамического воздействия на газонасыщенный 
породный массив. 


Техническая постановка задачи 


Пневмогидродинамическая обработка подземных нефтегазоносных массивов 
осуществляется путем высоконапорного нагнетания аэрированной жидкости через 
скважины, пробуренные с поверхности на продуктивный горизонт для изменения 
его состояния с целью интенсификации выхода полезного продукта. Кроме того, 
воздействие на уг ольный пласт осуществляется с целью его гидрорасчленения как 
метод предварительной дегазации, вызывающий активизацию выделения метана. 

На рис. 1 приведено схематическое изображение технологии воздействия на 


пласт через скважину, пробуренную с поверхности. 


$=0 
Ро 


- 


воздух 


ЖИДКОСТЬ 


Ва. 


Рисунок 1 — Схема пневмогидродинамического воздействия на пласт 
через скважину, пробуренную с поверхности 


Анализ результатов выполненных теоретических разработок 
в области математического моделирования процесса 


Сходство основных закономерностей процессов движения жидкостей и газов в 
угле и горных породах позволяет использовать для этих исследований разработки в 
области теории фильтрации в нефтегазовых коллекторах и грунтах. 

Основные модели напорной фильтрации жидкостей в горных породах описы- 
вают жесткий, упругий и нелинейно-упругий режимы для пористых и трещиновато- 
пористых сред. Жесткий режим фильтрации представляет собой движение несжи- 
маемой жидкости в недеформируемой пористой среде. В общем случае, при наличии 
в системе источников или стоков этот режим описывается уравнением Пуассона 
относительно давления [1]: 

ДР, (1) 


где Е} — давление жидкости; 
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Г — функция источника (стока), учитывающая свойства среды и фильтрующейся 
жидкости. 

Упругий режим соответствует нестационарному изменению давления слабо- 
сжимаемой жидкости в малодеформируемой пористой среде, когда жидкость и 
фильтрующая среда подчиняются закону Гука [1], [2]: 

2 ДАР. ©) 
01 
К 


п, в + В.) 
К — коэффициент проницаемости среды; 
п, — эффективная пористость; 


тд: = — коэффициент пьезопроводности; 


д - вязкость жидкости; 

В„.В. — коэффициенты сжимаемости соответственно жидкости и пористой 
среды. 

Уравнение (1) не учитывает изменение фильтрационных характеристик пори- 
стой среды в процессе ее деформации. В работе [3] на основе экспериментальных 
данных получена линейная зависимость проницаемости и пористости среды, плот- 
ности и вязкости жидкости от давления получено уравнение нелинейно-упругой 
фильтрации жидкости в пористой среде, по форме сходное с уравнением (2): 


= = хат [1+ о(Р — Р,)етааР}, (3) 


Ре 
где м = В, +В», — В,, 
В,.В,.В»„ - коэффициенты, характеризующие зависимость проницаемости 


среды и вязкости жидкости от давления; 
Р о - Пластовое давление. 


В каждом конкретном случае для выбора одного из уравнений (1) - (3) не- 
обходимо исследование степени влияния давления на характеристики пористой 
среды и жидкости и специфики фильтрационного потока. 

В реальных горных породах, представляющих интерес с точки зрения теории 
фильтрации, эффективная пористость представлена трещинами и порами с очень 
широким диапазоном зияния: от десятков ангстрем до миллиметров. 

В работе [2] предложена система уравнений, описывающая процесс движения 
жидкости в трещиновато-пористой среде: 


О" р) |, 

(ри) -9=0: | дР 
АЛЕ Ир) -В) 

О("р) и (1 

о о й ЭР (4) 
"2 ДР, = (В. +п,р) 2 —ок(Р, -Р 

4 = ра(Р, -Р), и 2 (Ве пр) д: ( 2 ‚) 


где 2... В.., В — коэффициенты сжимаемости среды 1, 2 и жидкости; 


индекс «1» относится к трещинам, индекс «2» —к блокам; 
9 — функция, характеризующая интенсивность обмена жидкостью между тре- 


щинами и блоками; 
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а — коэффициент, количественно характеризующий интенсивность обмена; 
р - плотность жидкости; 


у — скорость движения жидкости. 

Как и при выборе режима фильтрации, выбор пористой или трещиновато-по- 
ристой модели осуществляется на основе исследования соответствующих свойств 
горных пород и специфики решаемых задач. 

Кратко рассмотренные модели являются основными для исследования напор- 
ной фильтрации жидкости в горных породах. Существующие в настоящее время 
способы математического описания движения жидкости в угольных пластах основа- 
ны на этих моделях. 

Построение адекватных моделей на основе известных моделей подобных или 
сходных процессов в первую очередь предполагает выделение существенных специ- 
фических особенностей рассматриваемой задачи. 

К применению трещиновато-пористой модели для угольного массива также нужно 
подходить с учетом его специфических особенностей. В нефтегазовых коллекторах 
основной запас нефти (газа) сосредоточен в пористых блоках, кроме того, блоки имеют 
значительные размеры [2]. 

С точки зрения гидродинамики, особенностью процесса нагнетания жидкости в 
угольный пласт является то, что фильтрация происходит в среде, первоначально не 
заполненной жидкостью, тогда как большинство моделей напорной фильтрации по- 
строены для насыщенных жидкостью пористых сред. Учет этого факта приводит, во- 
первых, к нестационарности процесса фильтрации, во-вторых, к необходимости введе- 
ния дополнительных условий для вычисления координат фронта движущейся жидкости, 
что существенно усложняет решение и ограничивает область применения модели. 


Разработка математической модели 
пневмогидродинамического воздействия 


Уравнения фильтрации жидкости и газа имеют вид [5]: 


ОР оо 


0<х<=< ИВ: 5 
1 дх? бе х (г) (5) 
РР 
: В ра (6) 
Ох / 
| КР 
Е о | (7) 
| П к д „-_@ В.Т 
| °Е, | (чар) 


Здесь Р,Р, — соответственно давление жидкости и газа; 
ЕЕ, — модули сжимаемости жидкости и пористой среды 
п — абсолютная пористость; 

а, 6 — константы изотермы Лэнгмюра; 

К, — газовая постоянная; 


Т - абсолютная температура; 
КЕ) — координата границы раздела «газ — жидкость», определяемая из уравнения 


и АО. (8) 
Ш т бх 
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Использование этой модели позволило показать, что для расчетов движения 
жидкости в массиве с учетом его свойств с достаточной для практики точностью можно 
использовать формулы жесткого режима и получить методику определения пара- 
метров нагнетания. 

Постановка задачи для трещиновато-пористой среды, представляющая собой 
систему двух уравнений: 


КОР 0 

ЕЕ (9) 

д дх (61 

д 

—— =о@Р(И’ - И’). 10 

>. (0% ) (10) 

с условиями 

Р(0,)=Р.; (11) 
Р(х > 116) = 0: (х> В =0; 
(х ) (х ) (2) 
=;  К0)=0; 

Я __ Е ЭР 

Ш ип, Ох | _›, (13) 


где И’ — приращение влагонасыщения; 
Й’, — максимальный прирост влаги; 


а — коэффициент, характеризующий проницаемость блочной системы. 

Условия (11) — (13), хотя и записаны по-разному, имеют одинаковый фи- 
зический смысл, то есть выражают, в соответствии с законом Дарси, прямо про- 
порциональную зависимость скорости движения свободной поверхности потока 
жидкости от градиента давления на этой поверхности и позволяют явно вычислять 
ее координаты. Однако наличие в пласте зон с различной проницаемостью приводит 
к сложной, неправильной форме поверхности фронта, характеризующейся образо- 
ванием «языков» влажности. Пример моделирования процесса в данных условиях 
приведен на рис. 2, где показаны линии равного давления в условиях неравномерной 
проницаемости (в процентах от давления на скважине). 

100 


< 


30 
20 
10 


Рисунок 2 — Результаты моделирования воздействия через одиночную скважину, 
пробуренную с поверхности 


Необходимость явного вычисления координат свободной поверхности потока 
жидкости в условиях анизотропного пласта и взаимодействия встречных потоков чрез- 
вычайно усложняет, а иногда и делает практически неосуществимым решение задачи в 
плоскостной И пространственной постановке. К. настоящему времени не разработаны 
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модели, позволяющие автоматически учитывать границу области фильтрации без яв- 
ного ее выделения, дающие возможность эффективно решать задачи движения жидкости 
в сложном анизотропном пласте. 

В [5] рассматриваются процессы гидродинамического воздействия на угольные 
пласты как основа способа добычи угля и предотвращения газодинамических 
явлений. Здесь рассмотрена краевая задача для уравнения пьезопроводности: 


ОР), 9 Р(х,6) 


14 
[6] дх? р 


Задача решается на конечных промежутках времени Т, < <Т,, #1 = 0,.1,2...., и 


рассматривается процесс сброса давления. Аналитически получено решение в виде 
бесконечного ряда: 


р 2. а Е 
Рив = |+ + РР и 
0 =У = (Р.+Р. и) , 


7 с 7 
(15) 
с, . АлЛХ а 
у —м-— А. (=) —Р.. 
р-т”, + +-Р 


Полученное выражение позволяет выполнять расчет точечных значений В (х 5Ё ) С 
помощью компьютера. 


Выводы 


Установлены требования и предложена структура математической модели, 
учитывающей специфические особенности процессов, происходящих при внедрении 
газожидкостных смесей в трещиновато-пористую сплошную среду, в первую очередь — 
конечность скорости движения фронта рабочего агента. Модель предназначена к 
включению в состав системы автоматического управления процессом пневмогидро- 
динамического воздействия на подземные массивы. 

Показано, что необходимость явного вычисления координат свободной поверх- 
ности потока жидкости в условиях анизотропного пласта и взаимодействия встречных 
потоков чрезвычайно усложняет решение задачи в плоскостной и пространственной 
постановке, а возможным способом упрощения задачи является применение элемен- 
тов нечеткой логики. 
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Тре Майетайса! Моаейие т Фе Зуяет оГАиотайс 
Сопй’о/ о} Ргосез5 оГРиеитойуагодупатес АсНоп 


оп Оиаегетоипа Маубтле 


ТБе ргоет оЁ сопзёгасйоп оЁШе зузет о амютайс сопёго] о ргосез$ оЁ рпеитайс 
ап Нудгодупапис асйоп оп газ сощаше ипдегогойп соа|-госК таззуе аз а \’ау оЁ 
теме свапотпо оЁ 15 зваваз 1$ сопз14егеа. Еог ЧезстЬ ше, оЁ ргосез$ 1е тафетайс тодЕе|, 
Базе оп Фе зует ог попПпеаг едиайоп$ ул 1есбпоозтсаПу дейпеа Боипдагу соп@ю0п5$ 
1$ ргорозед. 

Те сБогсе оЁ рогоиз ап4 сгаск-рогоиз$ тоде| оР таззтуе таКез оп Фе Базе оЁ шуезй- 
саНоп оРсотезропате ргорег@ез оРпше госК ап4 зрестНс Геабалгез оЁ о1уеп ргоЫетпз. 

Тре сопзёгасНоп оЁ уечйе4 тоде!5 из Кпоууп то4е[$ оЁ зпаЙаг ог ПКепезз ргосеззез 
а Нт$Е Чите гедиге$ ЧзИполи$ те зрес!Йс рагасШагу оЁсопя4етпе адиезвоп. 

Егот Фе рошЕ оЁ Будгодупапс$ Ше рагасшагиу оЁ Фе ргосез$ оЁ Паша ритрие ш 
соа| зеат 15 фе соп4 оп, фа Ноулис 15 готе ш змтоипате, сощаиа$ по Паша БеЮге, Би 
1105 оЁпо4е]5 орритре ууаз сопзиасеа Гог Паш сощашед рогоцз тазз1уе. 

ТЫ$ Расе |еа45 5 ю поп Зап4ше о ргосез$ ап4 © песеззИу оЁ шсшАте ад опа! 
сопопз Гог са1сшайоп оЁ соог4таез оЁ Боипдапез оР Йоулпе Паша, Фа таКез зо[айоп 
тоге аси. 

ТБе ех15Нпо оР 7опез ог Чегет репегайпо ргоу14ез сотрИсае4 сигуед Юпп оЁ 
зигГгасе оЁ Паша НЙоулие Боипаахгу. 

Те сопзбасНоп оР тафетайс тоде| ум зрес!с Геашгез 15 песеззагу Гог шуезНоайоп 
ОР ргосеззез оЁ Йоулие Паша ап газ ш рогой$-сгаск таззуе ш соп4юп$ оЁ пщегасНоп оЁ 
тоушо веаплз. 

ТБе ппрокапё гое оЁР тафетайса| тодеПпс {фаКез расе Гог зооп оЁ фе рго Мет оЁ 
Впа! заое оЁ асйоп, рго\14ез зисВ аеогт, Фаё Верз ю тефисе пегайуе шЙиепсе юг 
есо]ор1са| соп@ опз. 


Статья поступила в редакцию 13.11.2012. 
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